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We investigate pseudogap phenomena realized in the BCS pairing model with a long but finite 
interaction range. We calculate the single-particle self-energy in all orders exactly in the temperature 
£SJ \ range where the superconducting fluctuation propagator is Gaussian-like. It is found that vertex 

» I . corrections to the self-energy, which are discarded in the previous studies, are crucially important 

for the pseudogap of the single-particle density of states in higher order calculations. 
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It is known that a pseudogap of the single-particle density of states (DOS) is a universal phenomenon observed in 
electron systems interacting with bosonic critical fluctuations.Eftj The interaction with a Gaussian fluctuation which 
arises near the instability toward a long-range order like antiferromagnetism, charge density wave, or superconductiv- 
ity results in a much enhanced damping of single-particle excitations, leading to a decrease of the DOS at the Fermi 
£j ■ level. Recently, it has been proposed by some authors that the pseudogap phenomena observed in high T c cuprates 
may be attributed to the superconducting nuctuation.QcJ They have succeeded in-e^hjjning .the pseudogap observed 
(*S in ARPES, NMR, and tunnel spectroscopy experiments. In the previous studies Jd'0'l3Ooil3'E3Efl the single particle 
self-energy is calculated by one-loop approximation (t-matrix approximation, diagrammatically expressed as Fig. 1(a)) 
or self-consistent t-matrix approximation (expressed as Fig. 1(b)). Although these calculations show qualitative agree- 
ment with experimental observations, there are some important differences between the results obtained by different 
approximation methods. For example, according to refs. 16 and 18, in the self-consistent one-loop approximation, the 
pseudogap of the DOS is destroyed near the Fermi level, and the quasiparticla.peak is restored. On the other hand, 
'-^J | such a restoration is not seen in the lowest-order one-loop approximationJ13'Ej Thus, it is desirable to examine the 
£5 ■ effects of vertex corrections, which are expressed as Fig. 1(c), and previously discarded without any justification of the 
approximation. Moreover, the previous studies suggest that the substantial decrease of the DOS at the Fermi level 
appears for the interaction stronger than the intermediate coupling. For this strong interaction, the validity of the 
t-matrix approximation is obscure. Motivated by these considerations, here, we investigate a simple model for which 
fS| , we can calculate all-order self-energy corrections in the temperature range where the superconducting fluctuation is 
k-* Gaussian-like. 

o 

oo 
m 



(a) * (b) 





(c) 



O 



O 
o 





XX7 



FIG. 1. Diagrams for the irreducible single-particle self-energy. The solid line represents the single electron Green's function. 
The wavy line represents the propagator of the superconducting fluctuation, (a) Lowest order one loop correction, (b) 
Self-consistent one loop approximation, (c) Some typical diagrams of vertex corrections which are not included in the diagrams 
(a) and (b). 
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The model Hamiltonian is given by 



H = J2 E k c\ a c ka - 1 V(q)B^q)B(q), (1) 

k,a q 

n^ni-^L (2) 

where B(q) — ^ k (^kCkic—k+q] , and d is the spatial dimensionH The first term of eq.(H) is the kinetic energy with 
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the energy dispersion Ek- The second term is the pairing attractive interaction with the interaction range l c . B(0) 
is a local annihilation operator of a Cooper pair with a structure factor For example, Cfc = 1 for s-wave pairing, 
£k = sin k x , sin k y for p-wave pairing, and £u = cosk x — cos k y , cos k x cos k y for cZ-wave pairing, and so forth. In the 
limit of l c — > +(}&. .the second term of eq. (|lj) is reduced to VB^(0)B(0). Then eq.(|l]) is the pairing model, which is 
exactly solvable.EacJ In this limit, the BCS mean field solution is exact, and the effects of fluctuation are completely 
suppressed. For finite but sufficiently large l c , the superconducting fluctuation is restored, and can be calculated 
systematically using the expansion in terms of We should bear in mind that although in this paper l£ is taken 
as a large parameter, Vl% = U must be a finite constant even in the limit of l c — ► +00 in order to keep the mean field 
transition temperature finite, and U is not a quantity of order 0(Z~). 

Superconducting fluctuation propagator p-|We first discuss the propagator of the fluctuating bosonic field. For this 
aim, following Krishnamurthy and ShastryJHj we introduce the auxiliary boson field <fi(q, r), and rewrite the interaction 
term of the action as 



iE Vj\F Jo iq ' T)B{<1, T) + h ' C] - (3) 



Integrating the fermion fields, we obtain the action for the superconducting order parameter fluctuation. 

S = E(* + ^ + ^ml^m^m 



+ C E $t>+q,m 2 +m s <i>k-q,m 1 -m 3 <i>k,m 1 <f>k', m2 +0{{l/lf) 2 ). (4) 
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Here, 4> q , m is the Fourier transform of <j>(q, r). t = (T-T c )/T c , a = 3C(2)7V(0)J7/(47rT c ), b = 7C(3)7V(0)C/^/(327r 2 r c 2 ), 
and c = 7C(3)iV(0)J7/(167r 2 ^T c 2 ) with T c the mean field transition temperature, and N(0) the unrenormalized DOS 
at the Fermi level. In the derivation of the quadratic term, we assumed that the renormalization of the chemical 
potential caused by the superconducting fluctuation is negligibly small. In a two dimensional (2D) case, at sufficiently 
low temperatures, the chemical potential may be strongly renormalized to-, a negative value, and eventually, the 
crossover to the Bose-Einstein condensation occurs at the zero temperature.Ej In this paper, we do not consider such 
a low-temperature region in the 2D case. We also ignore the momentum and frequency dependencies of the coefficient 
of the quartic term c. The momentum dependence is negligible, since it is of a higher order in The energy 

dependence is also negligible in the low energy region oj <C Ep. Note that c is of the order ~ Thus, in the limit 

of l c — * +00, the Gaussian term (the first term of eq.(|J)) dominates, and the propagator of the boson field is given by 



f 



bq 2 + t 



(5) 



This result was previously obtained by Krishnamurthy and Shastry@ In this limit, the BCS mean field solution is 
exact, and as a result, the single-particle self-energy correction vanishes. 

To compute the effects of superconducting fluctuation on the self-energy, we consider finite but very large l c . Then, 
we need to take into account the effects of the quartic and higher order terms of eq.(||). In the following, we consider 
only the quartic term, since it gives leading corrections. Applying the standard renormalization group method,c3 we 
obtain the one-loop scaling equations, 

^= 2 *+T JL T> (6) 



dl 1 + t 

dc b c 2 
di = {2 - d)c -—t) 



(2-^-77^- (7) 



Here, ao and bo are positive constants. From eq. (]?]), we see that for d > 2, the quartic term is irrelevant or marginally 
irrelevant. The main effect of the quartic term is the renormalization of the transition temperature, which is estimated 
from cq.(Q). The coefficients a and b of eq.(|^) are renormalized by higher order corrections than two loop. Thus, if 
the Ginzburg criterion for the smallness of fluctuations is satisfied by the initial values of the parameters, t, a, 6, and 
c, the Gaussian- type propagator eq.(|jJ) with the renormalized t, a, and b can be regarded as the exact propagator of 



2 



the superconducting fluctuation. Now we consider the case of t ^> bl c 2 in which, as will be shown later, the all-order 
self-energy corrections can be calculated exactly. In this case, the Ginzburg criterion is given by 



4cT e 



< t. 



(8) 



Since c ~ 0{l/lf), the condition eq.(||) is satisfied when l c is sufficiently large to satisfy t bl~ 2 . In the following, 
we consider only the situation wherein the unrenormalizcd coefficients a, b, c, and t satisfy the Ginzburg criterion 
eq. (j^) , and that eq. ([s]) with the renormalized parameters is a good approximation for the correlation function of the 
superconducting fluctuation. 
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FIG. 2. (a) Some typical self-energy diagrams in which all initial vertices are on the left of all final vertices. The symbol i 
represents the initial vertex, and / the final vertex, (b) Some examples of diagrams which are not included in £^ (e). 

All-order self-energy corrections - We expand the irreducible single-particle self-energy in terms of 



E fc (e) = ^4" ) (e). 



(9) 



Here, £^(e) is the n-th order term in 1//^. El (e) consists of two parts: one part is the contribution from the 
diagrams in which all initial vertices are on the left of all final vertices as shown in Fig. 2 (a), and the other part is the 

contribution from other irreducible diagrams as shown in Fig.2(b). We express the former one as Hi™"' (e). Carrying 
out the frequency summation analytically, we obtain 



{p} 



9k{q m ) 



■ ie + (-l)" 1 " 1 ^, 

qi,. ..,q n m=l K ~ / ,< =1 H 



i E 7 =i C '(9j) s g ne 



n 



iE + (-l)m'-iE , + iE™ i C(q Pj ) S gne 



(10) 



where gk(q) — T£lV(q)/(t + bq 2 ), and C(q) = (t + bq 2 )/a. {P} = {P 1 ,P 2 , ■~,P n } is generated by a permutation of 
{1,2, ...,n} under the rule that permutations between numbers with different parities are not allowed. For example, 
{P} = {1,2, 3}, {3, 2,1} for n = 3, and {P} = {1, 2, 3, 4}, {3, 2, 1, 4}, {1, 4, 3, 2}, {3, 4, 1, 2} for n = 4, and so forth. 
J2{p} 1S trie summation over these permutations. 
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In general, we cannot carry out g-summation of cq. (|l0|) analytically. However, for t b/l%, eq.([To|) can be greatly 
simplined. As we will see later, this case is very important for us. In this limit, we can replace g-summation of 
eq.(l(]), jf J2 q i w ith (Vc) E 9 ^9,o> since the momentum transferred is restricted within q < l/l c <C kp, for which 

gk{q) ~ <7/c(0), C(q) ~ C'(O), and E^- q ~ Ek- In this approximation, all crossing diagrams of e£™^ (e) give the same 
contribution as its non-crossing diagrams. This means that vertex corrections which are not included in self-consistent 
t-matrix approximation are not negligible. The total number of these diagrams, N(n), is m!(m + 1)! for n = 2m + 1, 
and (ml) 2 for n — 2m. Thus, we can rewrite cq. (|l0|) as, 

S fe (e) = N(n)gl i£ + (_y n _ lEh + inC (0) sgas li (ie + (-l)™-i£ fc + imC(O)sgne) 2 ' (H) 

where g k = gk(0)/lf. This expression has the same fiprm as that obtained by Sadovskii for one-dimensional fermion 
systems interacting with static Gaussian fluctuation.^ Then, we can apply Elyutin-Sadovskii's combinatorics method 

to evaluate all diagrams. According to Sadovskii, the contributions from T,^ (e) — sl n ' (e) are expressed in terms of 

^i" 1 ^ ( e ) with m < 7i.@ Using Elyutin-Sadovskii's method, we can write the self-energy in all orders as a continued 
fraction form, 

s fc(e) = To\~ > ( 12 ) 



Si(k,e) 



S 2 (k,e) V 



S 3 (k,s) 



S 4 (k,s) 



where S m (k, e) = ie + (— l) m_1 i?fc + irnC(0)sgne, and v(m) — [(m+ l)/2], [...] is Gauss's symbol. From this analytical 
expression, we obtain the single-particle DOS, p(e) = — 1/ir J2k I m [l /( £ — — Ef (e))] . We calculate p(e) numerically 
in the case of s-wave pairing = 1, replacing ^ fc with -/V(O) J de k - The result is shown in Fig. 3(a). The parameters 
are chosen so as to satisfy the condition eq.(^) and bl~ 2 <C t. We see that the pseudogap appears near the Fermi level. 
We also show the energy dependence of the self-energy in Figs. 3(b) and 3(c). It is seen that the scatterig rate, the 
imaginary part of the self-energy, is much enhanced near the Fermi JLea^el, and that the real part has a positive slope 
for e ~ 0. These features are characteristic of the pseudogap stateBtM A similar result is obtained for other pairing 
symmetries, such as p-wave, <i-wave, and so forth. 
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FIG. 3. (a) The single particle DOS p(e) plotted as a function of energy. The vertical axis is renormalized by the bare DOS 
po- g = 1, C(0) = 0.1 for the solid line, g = 1, (7(0) = 0.5 for the dotted line, (b) The imaginary part of the self-energy 
at the Fermi momentum plotted as a function of energy. The same parameters as Fig. 3(a) are used, (c) The real part of the 
self-energy at the Fermi momentum plotted as a function of energy. The same parameters as above are used. 



Analytical results for single-particle damping - The pseudogap behavior is due to the large enhancement of the 
damping, which is caused by strong scattering with superconducting fluctuation. Here, we will see how this enhance- 
ment is affected by the vertex corrections to the self-energy which were discarded in the previous studies. For this 
purpose, we calculate analytically the damping of electrons at the Fermi level, 



T.nE« (0) = ( 

cr(0) + v{2)9k 



2(7(0) + * 



3(7(0)' 



4(7(0) + 

The limiting value of this continued fraction is obtained from the recursion relation, P n +\ = (n + 1)(7(0)P„ + v(jiJf. 
2)flP n -i, Qn+i = (n + l)(7(0)Q n + v(n + 2)pQ„_ 1 , P = 0, Q = 1, P-i = 1, Q-i = 0, and 7 = lim^,. P n /Q n H 
Up to a numerical constant, we have 7 ~ gu/C(Q) = aTU/(TT d l^t 2 ). For y/aU/Jn^lf) 3> f, the damping of electrons 
is much enhanced, ~ T/t 2 3> T. Since we consider the case of t 3> &/'ci for which the Ginzburg criterion eq.(||) is 
also satisfied, the temperature range where the above result is applicable is b/l 2 <C t <C y at/ /(Tr d l d ). It is noted that 
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the result of all-order calculation gives the same asymptotic behavior as that obtained by the lowest-order one-loop 
approximation which is also given by ~ g/C(0) ~ T/t 2 . This implies that higher order corrections cancel each other, 
and the lowest-order term gives the leading contribution. This cancellation is not due to the presence of the Ward 
identity. We do not have the Ward identity which relates a self-energy to a vertex in particle-particle channel. Actually, 
the cancellation takes place among different order diagrams. The total number of the diagrams which give dominant 
contribution in the n-th order increases as 0((n\) 2 ). As a result, higher order contributions cancel with lower order 
contributions effectively. We compare this result with that obtained from self-consistent i-matrix a ppr oximation. The 



self-energy in the self-consistent approximation is derived by substituting v(m) — 1 for all m of eq.(12). The damping 
at the Fermi level, 7 SC , is given by 



c(o) + Mm ) 

2C(0) + ^ 



3C(0) 1 ''''' 



4(7(0) + 

where J n (x) is the n-th Bessel function. For sufficiently small t (large y/gk/C(0)), eq.([T^) is well approximated 
by \/gk = KfcC^T 1 /(^c 7r<i )] 1 ^ 2 ~ 1/V*' This is much suppressed compared with the result of all-order calculations. 
For t > b/l 2 , <C [t^UT/^blc -2 )] 1 / 2 . Thus, in the case of d — 2 or 3, the damping 7 SC cannot be enhanced 
sufficiently under this condition. This result presents a clear contrast with the all-order result. The self-consistent 
i-matrix approximation, which neglects vertex corrections underestimates crucially the effects of strong scattering by 
critical fluctuations, and may cause a fatal error not only quantitatively but also qualitatively in a pseudogap state. 
These results instructively show that the lowest-order one-loop calculation is more reliable than the self-consistent 
one-loop treatment, and furthermore imply that the lowest order one- loop treatment of the self-energy corrections 
may be a fairly good approximation because of the cancellation between higher order corrections. 

In summary, we have studied the pseudogap of the single particle DOS in the BCS pairing model with a long but 
finite interaction range using a well controlled approximation method. We have derived the all-order single-particle 
self-energy exactly in the temperature range where the superconducting fluctuation propagator is a Gaussian form. It 
has been revealed that the vertex corrections to the self-energy are essentially important in higher order calculations. 
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supported by a Grant-in- Aid from the Ministry of Education, Culture, Science, and Technology, Japan. 
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